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Injection or optical generation of electrons and holes in graphene bilayers (GBLs) can result in the 
interband population inversion enabling the terahertz (THz) radiation lasing. The intraband radia¬ 
tive processes compete with the interband transitions. We demonstrate that remote doping enhances 
the indirect interband generation of photons in the proposed GBL heterostructures. Therefore such 
remote doping helps surpassing the intraband (Drude) absorption and results in large absolute val¬ 
ues of the negative dynamic THz conductivity in a wide range of frequencies at elevated (including 
room) temperatures. The remotely doped GBL heterostructure THz lasers are expected to achieve 
higher THz gain compared to previously proposed GBL-based THz lasers. 
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I. INTRODUCTION 

The population inversion created by injection or opti¬ 
cal pumping in graphene layers (GLs) and graphene bi¬ 
layers (GBLs) [lj results in the negative dynamic conduc¬ 
tivity in the terahertz (THz) range of frequencies @, Jj 
and enables the graphene-based THz lasers HH. Dif¬ 
ferent research groups have demonstrated the THz gain 
in the pumped GLs (see, also Ref. 0 )- GLs 

and GBLs can serve as the active regions of the THz 
lasers with the Fabry-Perot resonators, dielectric waveg¬ 
uides, slot-lines including the plasmonic lasers [§!■ fl9l - [24l ] . 
The GL and GBL heterostructure lasers can operate in a 
wide frequency range, including the 6 to 10 THz range, 
where using materials like A 3 B 5 is hindered by the op¬ 
tical phonon effects (25l - f27j . The THz lasing in GL and 
GBL structures is possible if the contribution of the in¬ 
terband radiative transitions to the real part of the dy¬ 
namic conductivity Re a in a certain range of frequen¬ 
cies surpasses the contribution of the intraband radia¬ 
tive processes associated with the Drude losses. The in¬ 
terband transitions include the direct (with the conser¬ 
vation of the electron momenta) and the indirect tran¬ 
sitions (accompanied with the variation of the electron 
momentum due to scattering). The limitations imposed 
by the momentum and energy conservation laws allow 
only for the indirect intraband radiative transitions. In 
the GLs and GBLs with sufficiently long carrier momen¬ 
tum relaxation time r, the direct interband radiative 
transitions dominate over the indirect intraband tran¬ 
sitions and Re a < 0. Recently 0,0, we demonstrated 
that in GLs and GBLs with primarily long-range disorder 
scattering, the indirect interband transitions can prevail 
over the indirect intraband transitions leading to fairly 


large absolute values of the negative conductivity | Reoj. 
These values could exceed the fundamental limits for the 
direct transitions, which are <tq = e 2 /4 h for GLs and 
2 oq for GBLs [l[ (here e is the electron charge and h is 
the Planck constant). As a result, in the GL and GBL 
structures with a long-range disorder one can expect that 
the condition Re a < 0 could be fulfilled in a wider range 
of the THz frequencies, including the low boundary of 
the THz range. When the momentum relaxation is asso¬ 
ciated with long-range scattering mechanism, moderate 
values of r might be sufficient for achieving Re a < 0, es¬ 
pecially in GBLs where the density of states (DoS) p(e) as 
a function of energy e near the band edges (p(e) — const ) 
in the latter is considerably larger than in GLs (where 
p(e) oc |e|). 

This paper deals with the GBL heterostructures with 
a remote impurity layer (RIL) located at some distance 
from the GBL plane and incorporating donors and ac¬ 
ceptors or donor and acceptor clusters. We consider the 
GBL-RIL heterostructure as an active region of an in¬ 
jection THz laser. We show that the electron and hole 
injection from the side n- and p-contact regions into the 
GBL-RIL heterostructures leads to the population inver¬ 
sion enabling the enhanced negative dynamic conduc¬ 
tivity, which, in turn, can result in lasing in the THz 
range. Figure 1 shows the schematic view of an injec¬ 
tion laser based on the GBL-RIL heterostructure at the 
reverse bias voltage V. The remote (selective) doping 
can markedly affect the scattering of electrons and holes 
(see, for example 0,0) and, hence, the indirect radia¬ 
tive transitions. The GL and GBL heterostructures using 
remote doping were already fabricated and used for de¬ 
vice applications 0, [0. The heterostructures with the 
selectively doped GBL and RIL demonstrate the follow- 
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FIG. 1: Schematic view of an injection laser based on GBL- 
RIL heterostructure with the GBL doped with acceptors and 
the RIL doped with both donors and acceptors and with the 
chemically doped n- and p- injectors. 


ing distinct features. First, the net electron and hole 
densities determined by both the pumping and RIL and 
the electron and hole quasi-Fermi energies, ep.e and £f.Id 
are generally not equal. Second, relatively large scale 
of the donor and acceptor density fluctuations result in 
smooth (long-range) spatial variations of the potential in 
the GBL created by the RIL, separated by distance d 
from the GBL. This could result in the electron and hole 
scattering with a relatively small momentum change and, 
consequently, increase in the contribution of the indirect 
interband transitions in comparison with the interaband 
transitions. Thus, the main role of the RIL is not to 
induce extra carrier in the GBL (as in high-electron mo¬ 
bility transistors) but to provide the specific mechanism 
of the electron and hole scattering reinforcing the indi¬ 
rect interband radiative transitions. We calculate the 
spectral characteristics of real part of the dynamic con¬ 
ductivity Re a in the pumped GBL-RIL heterostructures 
as a functions of the doping level and the spacer thickness 
d and compare the obtained characteristics with those of 
the undoped heterostructures. This analysis reveals the 
conditions for the selective doping enabling the negative 
dynamic conductivity in a wide frequency range, includ¬ 
ing the frequencies of a few THz. This might open new 
prospects of GBL heterostructures for the efficient THz 
lasers based on the GBL heterostructures. 


II. DEVICE MODEL 

We study the GBL-RIL heterostructures pumped via 
the injection from the side p- and n- contacts shown in 
Fig. 1. In the very clean pumped GBLs (or GLs) with a 
fairly high carrier mobility the intraband (Drude absorp¬ 
tion) at the radiation frequencies above one THz should 
not play a dominant role. Hence, adding the selective 
doping might not be needed. However, in the GBLs with 
a relatively low carrier mobility caused by unavoidable 
(residual) charge impurities and imperfections, the RIL 
induced scattering could result in a substantial compen¬ 
sation of the Drude absorption and, therefore, enhance¬ 
ment of the negative dynamic conductivity In the fol¬ 
lowing, for definiteness, we consider the GBL-RIL het¬ 
erostructures with the GBL doped with acceptors (with 
the density E^) and the RIL formed by partially com¬ 



FIG. 2: Energy band diagram of a pumped GBL (left panel) 
and energy dependency of its DoS (right panel). Arrows cor¬ 
respond to direct and indirect interband transitions with the 
photon emission (wavy arrow ” 1” and wavy + smooth arrows 
” 2”) as well as to indirect intraband transitions with the pho¬ 
ton absorption in the conduction and valence bands (wavy + 
smooth arrows ”3” and ”3*”, respectively). 


pensated donors (with the density E d,r) and acceptors 
(with the density T,a,r)- Hence E r>,r — E a,r — E^. 
We further assume that in these heterostructures, the 
RIL might comprise the clusterized acceptors and donors 
(with the correlated charge defects forming the charged 
clusters with the charge Z c e > e and the characteris¬ 
tic size l c , where e is the electron charge). In the lat¬ 
ter case, we assume that the cluster density is equal 
to (Sa,,r + T,r>,r)/Z c . Such GBL-RIL heterostructures 
are quasi-neutral even in the absence of the pumping. 
Pumping results in the formation in the GBL of the two- 
dimensional (2D) electron and hole gases with the equal 
densities E e = E^ = E and the quasi-Fermi energies 
£p, e = £F,h = £f- Figure 2 shows the GBL band di¬ 
agram under pumping. In Fig. 2, the arrows indicate 
different radiative transitions: vertical interband transi¬ 
tions, indirect intraband transitions accompanied by the 
electron scattering on the impurities, phonons, and holes 
(leading to the Drude absorption of radiation), and the 
interband indirect transitions. 

The transverse electric field arising in the spacer due 
to the features of the doping can lead to the local open¬ 
ing of the energy gap l], |34|. However, as shown be¬ 
low, for the doping levels under consideration, the en¬ 
ergy gap expected to be small, and the band open¬ 
ing and and the e nerg y spectrum nonparabolicity will 
be disregarded [3J-|36|. The energy spectrum is as¬ 
sumed to be e p = ±p 2 Vw/ 7i = ±R 2 / 2 m, 0 , where 

71 ~ 0.35 — 0.43 eV is the inter-GL overlap integral in 
GBLs, Vw — 10 8 cm/s is the characteristic velocity, and 
to = 7 i/ 2 ujy is the effective mass of electrons and holes. 
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III. MAIN EQUATIONS 


The real part of the net dynamic (THz) conductivity 
in the in-plane direction, Re a, of GLs and GBLs com¬ 
prises the contributions from the direct (vertical), Re er^, 
and two types of the indirect, Re cr™^ er and Re 
transitions, respectively: 


Rea = Re(a d + a™ t d er +a™ t d ra )- W 

The first term in Eq. (1) corresponds to the transitions of 
the type ” 1” in Fig.2, the second term corresponds to the 
type ” 2”, and the third term - to the types ” 3” and ” 3*” 
transitions, respectively. The last term in the right-hand 
side of Eq. (1) corresponds to the processes responsible 
the Drude absorption. 

Due to high frequencies of the inter-carrier scattering 
under sufficiently strong pumping, the energy distribu¬ 
tions of the pumped carriers are characterized by the 
Fermi distribution functions f v (e p ) and / c (e p ), where e p 
is the dispersion relationship for the 2D carriers in GBLs, 
with the quasi-Fermi energy, £p — Trh 2 Y,/2m and the ef¬ 
fective temperature T. 

The calculations of the GBL dynamic conductivity in 
the in-plain direction associated with the direct inter¬ 
band transitions accounting for the GBL energy spec¬ 
tra, following the approach developed in {33, HU (see, 
also 1 ), yield: 


Here g = 4 is the spin-valley degeneracy factor, £ p , 
v p = 2pvyy/"/i = p/m is the velocity of electrons and 
holes with the momentum p in GBLs, V (q) is the q-the 
Fourier component of the scattering potential, Mpp, = 
(1 + AA' cos 20 pp /)/2 is the overlap between the envelope 
wave functions in GBLs, and A = ±1 is the index indicat¬ 
ing the conduction and valence bands. The extra factor 
of 2 in Eq. (3) accounts for the intraband electron and 
hole contributions. 

Equation (3) yields the well-known result for the dy¬ 
namic Drude conductivity in the high-frequency limit 
r = H/t < hu> < £f, where F is the broadening of the 
carrier spectra. 


IV. SCATTERING MECHANISMS 


For the GBL-RIL hererostructures under considera¬ 
tion, we have the following formulas: 


|U(fiq )| 2 


27r e 2 l 2 
_n(q + q T F)_ 


[£a + Zrcl exp(—2gd)], (5) 


for the discrete acceptors in the GBL and the discrete 
acceptors and donors in the RIL and 


Re 


££= (fc. + 2 7l ) taiih /fa — 2 £f A 
c TQ (ftw + 7i) \ 4T J 


~ 2 tanh 


luo — 2ep 
4 T 


( 2 ) 


Applying the Fermi golden rule for the indirect intra- 
and interband electron radiative transitions, we obtain 
at the following formulas for the respective components 
of the dynamic conductivity: 


Re 


' ind 

a Q 


4 TTff 
hu 3 


Y I^(P-P')| 2u pp'( v p' - Vp) 2 6(/!W+£p -£p') 


p,p' 


1 + exp[(£p - £f)/T] 1 + exp[(£p/ - £ F )] 


, ( 3 ) 


Re 


' ind 


2 tt 5 

Hlo 3 


Y l^(P-p')| 2 Wp p /( v p' +v p ) 2 5(fiu;-£p -£ P 0 


P,P' 


1 + exp[(£ p - £ F )/T] 1 + exp[(£ p / + e F )/T} 


( 4 ) 


\V{hq)\ 2 


2ire 2 l 2 
n{q + 9 tf). 

x [Ex + Z c Y, RC l exp(- 2 qd 


q 2 l 2 /2)}, ( 6 ) 


for the point charged defects (acceptors) in the GBL and 
the correlated clusterized acceptors and donors in the 
RIL, respectively. Here q = |p-p'|, £_r cl = £ a,r+'Z D ,r 
is the net density of the charge impurities of both types 
located in the RIL (therefore the cluster density in the 
’’clusterized” RIL is equal to T, R cl/Z c with £ F cl > 
Ex), k is the effective dielectric constant (which is the 
half-sum of the dielectric constants of the materials sur¬ 
rounding the GBL, i.e., of the spacer and the substrate), 
qTF = (4 e 2 m/h 2 K) is the Thomas-Fermi screening wave 
number [38j-|40(, which is independent of the carrier den¬ 
sity. Deriving Eq. ( 6 ), we have taken into account that 
the density of the charged clusters in the clusterized RCL 
is Z c -times smaller than the donor and acceptor densi¬ 
ties, but the cross-section of the carrier scattering on the 
clusters is proportional to Z 2 e 4 . The factor exp(—2 qd) in 
Eqs. (5) and ( 6 ) is associated with the remote position of 
a portion of the charged scatterers (see, for example 0). 
Equation (7) corresponds to the Gaussian spatial distri¬ 
butions of the clusters that leads to the appearance of the 
factor exp(— q 2 l 2 /2) (compare with Refs. [291. 1411 - 14^ 1. 

One can see that the scattering matrix element explic¬ 
itly depends on the background dielectric constant k and 
via the dependence of on k. 
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FIG. 3: Frequency dependences of Re + aZT) cal¬ 

culated for GBL-RIL heterostructures with spacer thickness 
d = 15 nm, different dielectric constants (k = 10,15, and 20), 
and different doping parameters D\ (a) D = 0, (b) D = 10, 
and (c) D = 20. 



FIG. 4: The same as in Fig. 3 but for doping parameter 
D = 20, different spacer thicknesses (d = 10,15, and 25 nm), 
and different dielectric constants: (a) k = 10, (b) k = 15, and 
(c) k = 20. 


V. NET THZ CONDUCTIVITY 


In the following we consider the most practical case of 
the temperature, frequency range, and pumping in which 
H/t < hw,T < Sf■ For T = 300 K, the frequency range 
uj/2tt = 1 — 10 THz, and Sf > 50 meV, these inequalities 
are satisfied when r > 0.05 — 0.5 ps. 

For GBL structures without the RIL (J^ril = 0) 
and the totally screened charges in the GBL, in which 
Re aZd™ = and Re a^T = , invoking 

Eqs. (3) - ( 6 ), one obtains 


Re er 


intra 

ind 


" 2 
mr 



( 7 ) 


a result, we arrive at the following simplified formulas for 

Re<cr + Re ^ r ): 


t> _ ( „intra i winter \ 
Refund + a ind ) 


to _intra . 

- Ke a ind < 



dqq 


x [1 + Dexp(— 2qd — q 2 l 2 /2)] 
qq (Q 2 max -Q 2 min )(d/<iTF + 1) 2 
[1 + D exp(—2 qd - q 2 l 2 / 2)] ) 
{qhax ~ 9mi n ){q/qTF + l) 2 ) ' 


where u>r> = \/4 HY,/mT = y/8 £f/fHt is the Drude fre¬ 
quency. This quantity characterizes the relative strength 
of the Drude processes (absorption). 

As the above takes place, Eqs. (3) and (4) lead to 


Re < 7 ™*™ 4e F ' 


( 8 ) 


To avoid cumbersome calculations and obtain trans¬ 
parent formulas suitable for a qualitative analysis, we 
use the mean value theorem for the integrals over dq. As 


Here D = Z c {Y* R cl/^a) = Z C (Y, D ^ R + T,a,r)/{^d,r — 
Zja,r) — Z c {jd D ,R + X a ,r)/Xa > Z c > 1 is the dop¬ 
ing parameter (Z c = 1 and l = 0 for the RGL with 
the point charges (acceptors and donors) and Z c > 1 
and l = l c in the case of RIL with the clusterized ac¬ 
ceptors and donors), hQ m ax — 2 -y/ 2 m£ F (l — hui/Aep), 
ftQmin — 2- v /2to£ F (&u/4£ F ), Hq max ~ 2 y/mhuj,q min ~ 0. 
The latter quantities follow from the conservation laws 
for the indirect intraband and interband radiative tran¬ 
sitions. In particular, Eq. (9) gives rise to 
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FIG. 5: Areas in the D—d plane where Re (cr^ T ’ a +crJ"^ er ) > 0 
(areas below lines corresponding to given frequencies) and 
Re + o'ind 6 ’') < 0 (above these lines) at different values 

of dielectric constant k: (a) k = 15 and (b) k = 20. 


Re (cr 


intra 
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) 


Retr 


intra 

ind 


~ l 


r q max [l + Dexp(-2qd-q 2 l 2 /2 )] 

i L. " (,/«t F+ d* 

2 rQmax [1 + Dexp(— 2qd — q 2 l 2 /2)] 
Jo min " (q/qtf + 1) 2 


( 10 ) 


For sufficiently pure GBLs and relatively strong RIL 
doping (with nearly complete compensation of the donor 
and acceptors), the doping parameter parameters D can 
be fairly large. In this case, the contribution of the carrier 
scattering on the remote impurities to the net scattering 
rate can be dominant. In such a case, the absolute value 
of the net THz conductivity |Recr| can markedly exceed 
the contribution solely from the direct interband transi¬ 
tions, i.e., |R.ecr| > 2 ctq. 

Considering Eqs. (1), (2), and (9), we arrive at the 
following formula for the net THz conductivity: 


Re a 

2 CTQ 


-*(t) 


LUd 
uj 
fuj 
4 ef , 


rQn 


dqq 


[1 + D exp(—2 qd — q 2 l 2 / 2)] 


(Qmax ~ QmiJid/dTF + l ) 2 

( 11 ) 


Qmin \^max mm 

Qmax [1 + D exp(—2 qd — q 2 l 2 /2 )] 
mi „ qq (Qla X -V 2 m MqTF + l ) 2 


VI. RESULTS AND ANALYSIS 

Using Eqs. (9) and (10), one can analyze the condi¬ 
tions when the interband radiative processes surpass the 



Spacer thickness, nm 

FIG. 6: Threshold frequency u>o/2f corresponding to 

Re (o~i£d ra + <J™T) = 0 as a function of the spacer thick¬ 
ness d calculated for D = 20 and different values of dielectric 
constant k. AS in Fig. 5, areas above the pertinent line are 
related to Re (a^ ra + (Tind er ) < 0. 


intraband one’s, i.e., when Re (cr.”^ ro + cr-™^ er ) < 0. 

For the practical GBL-RIL heterostructures with the 
Hf 02 spacer and the substrate material like Si and SiC, 
the effective dielectric constant (which is the half-sum of 
the dielectric constants of the spacer and the substrate) is 
in the range k = 10 — 20. In the calculations we assume 
set to = 4 x 10 -29 g and Ef = 60 meV (E ~ 2.5 x 
10 12 cm -2 ). 

Figures 3 and 4 show the dependences of Re (cr™^ ra + 
a ind r ) on lh e frequency / = u}/2?r calculated using 
Eq. (9) with l = 0 (uncorrelated impurities in the RIL) 
for different values of the doping parameter D , the spacer 
thickness d, and effective dielectric constant k. As seen 
from Fig. 3(a), in the absence of the RCL (D = 0), the 
indirect intraband transitions prevail over the indirect 
interband transitions. In such a case, the indirect tran¬ 
sitions lead to the photon absorption in the entire fre¬ 
quency range, although the photon emission due to the 
interband transitions partially compensates the Drude 
absorption. This compensation is enhanced in the het¬ 
erostructures with the doped RIL, i.e., for D > 0 [see 
Figs. 3(b) and 3(c)]: an increase in the doping param¬ 
eter D leads to the appearance of the frequency range 
where Re (cr' L £ l t ( [ a + cr-^ er ') < 0. The latter range becomes 
wider (it starts from lower frequency) with increasing di¬ 
electric constant k. These features of the transformation 
of the Re (o--"^ ra + cr-"^ er ) frequency dependence are at¬ 
tributed to reinforcing of the long-range carrier scattering 
(with increasing D ) and smoothing of the scattering po¬ 
tential (due to a stronger screening at higher k). Both 
these effects result in an increase of the relative role of 
the indirect interband transitions. The scattering poten¬ 
tial smoothing associated with an increase in the spacer 
thickness d also beneficial for the indirect interband tran¬ 
sitions is seen in Fig. 4 as well. 

Figure 5 shows the D — d plane where Re (cd"))™ + 
a ind er ) > 0 (below the D versus d lines) and Re (a™*[ a + 
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Frequency, THz 


FIG. 7: Frequency dependence of the normalized net THz 
conductivity Re<r/2crQ of GBL-RIL heterostructures with k = 
15 and different values of doping parameter for: (a) d = 15 nm 
and cud/ 27 t = 2.49 THz (r = 1 ps), (b) d = 25 nm and 
= 2.49 THz(r = 1 ps) , and (c) d — 25 and cud/27t = 
3.51 THz (r = 0.5 ps). Lines marked by circles correspond 
to normalize net THz conductivity without RIL (D = 0). 
Dashed lines show normalize THz conductivity due to solely 
direct interband transitions Re<Td/2o-Q. 


a ind r ) < 0 (above these lines) calculated considering 
Eq. (10 )at different values of k and frequency / = ui/2n. 
As seen, the ranges of parameters D and d where the 
quantity Re (cr-"^™ + tT-”^ er ) is negative are markedly 
wider in the case of relatively large dielectric constant k 
and the frequency /. 

Figure 6 shows the threshold frequency ujq/2tt corre¬ 
sponding to Re (a"///[ a + <r”^ er ) = 0 as a function of the 
spacer thickness d calculated for D = 20 and different 
values of k. 

Figure 7 shows the frequency dependences of the nor¬ 
malized net THz conductivity Re a/2a q calculated us¬ 
ing Eq. (11) for T = 300 K and different values of 
the doping parameter D 1 dielectric constant k = 15, 
and Drude frequency £Id/ 27t. The Drude frequencies 
are chosen to be wd/27t = 2.49 and 3.51 THz. At 
£f = 60 meV, this corresponds to the carrier momen¬ 
tum relaxation time r = 1 ps and 0.5 ps, respectively. 
The latter values of r corresponds to the carrier mobil¬ 
ity in the GBL without the RIL, i.e., with D = 0, equal 
to fi = 40, 000 cm 2 /V s and 20,000 cm 2 /V s. The large 
values of D imply rather strong doping of the RIL. For 



Frequency, THz 


FIG. 8: Frequency dependence of the normalized net THz 
conductivity Re a/2aQ of GBL-RIL heterostructures with 
clusterized impurities, k = 15, cj/27t = 3.51 THz, d = 15 nm, 
D = 60, and different values of the cluster size l c . 


example, if T,a = 1 x 10 11 cm 2 (that is consistent with 
the above value of r), the parameter D is equal to 60 if 
Y, D , R = 3.05 x 10 12 cm " 2 and H D , R = 2.95 x 10 12 cm" 2 , 
g. As seen in Fig. 7, the net THz conductivity is neg¬ 
ative in a certain range of frequencies (above 1.5 THz) 
even for the case D = 0, i.e., even in the GBL structures 
without the RIL. This is because although the indirect 
intraband transitions (Drude transitions) prevail over the 
indirect interband transitions, the negative contribution 
of the direct interband transitions enable the negativity 
of the net conductivity. The latter is because the contri¬ 
bution of the direct transitions decreases with increasing 
frequency much slowly (see the dashed line in Fig.7) than 
the contribution due to the Drude processes. However, 
in the GBL-RIL heterostructures with sufficiently strong 
doping (large D) the absolute value of the net THz con¬ 
ductivity is pronouncedly larger than in the similar het¬ 
erostructures without the RIL. Moreover, the substan¬ 
tial contribution of the indirect interband transition can 
result in |Re cr| markedly exceeding its values |Re ad\ 
associated with solely direct interband transitions (see 
the curves s for D = 40, and 60 and the dashed lines in 
Fig. 7) and even exceed the maximum value of the latter 
2 ctq (see the curve for D = 60). The comparison of the 
frequency dependences in Fig. 7 for D = 0 and D = 60 
reveals that the RIL can provide a substantial increase 
in the absolute value of the THz conductivity [up to 4 - 
6 times increase in the range oj/2-k = 3 — 5 THz ,com¬ 
pare the curves for D = 60 and D = 0 in Fig. 7(c)] and, 
hence, the corresponding enhancement of the laser THz 
gain (for the parameters under consideration). Such an 
increase can be even more large at somewhat lower fre¬ 
quencies. 

The comparison of Figs. 7(b) and 7(c) also demon¬ 
strates that at larger Drude frequencies, i.e., at shorter 
carrier momentum relaxation time r (and the same pa¬ 
rameter D ), the sag of the net THz conductivity fre¬ 
quency dependences becomes deeper. This is attributed 
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to an increase of the intensity of the indirect processes 
(both intraband and interband) when r decreases. 

The plots in Fig. 8 were obtained using Eq. (11) with 
l = l c ^ 0. Figure 8 shows the frequency dependences 
of the normalized net THz conductivity in the GBL-RIL 
heterostructures with clusterized (correlated) impurities 
in the RIL. As seen from Fig. 8 , for relatively small 
spacer thickness (d = 5 nm), the clusterization leads to 
an additional enhancement of the indirect interband tran¬ 
sitions, reinforcing the negative THz conductivity, partic¬ 
ularly with an increasing cluster size. However, at rela¬ 
tively large spacer thicknesses (d > 10 nm), an increase 
in the absolute value of the net THz conductivity with 
increasing size of the clusters is insignificant, since in this 
case, D is determined not only by the doping of the RIL 
but the cluster charge Z c . This means that the chosen 
value D = 60 can correspond to different densities 
and S ^4 r depending on Z c . The latter, in turn, depends 
of the degree of the compensation of donors by acceptors 
in the clusters. 


VII. DISCUSSION 

If the RIL doping does not compensate the acceptor 
system in the GBL, the electron and hole quasi-Fermi en¬ 
ergies £F,e and £F,h are generally not equal to each other 
even in the pumping conditions. However, at sufficiently 
strong pumping when the density of the injected (or op¬ 
tically pumped) carriers £ > AE = T,f>,a — Z,d,a — ^d,a, 
Eqs. (2), (9), and (10) are still approximately valid until 
tuo/2 < min{£ Fie ,£f i?l }. 

The above model disregarded the opening of the band 
gap in the GBL under the transverse electric field aris¬ 
ing due to the impurity charges in the RIL and GBL 
Despite a fairly complex pattern of th e op ening gap in 
GBLs by the electric field and doping |4l|, [42| , the per¬ 
tinent energy gap A g can roughly be estimated in the 
manner as it was done in (35| . In the absence of the 
electron and hole injection from the side contacts, the 
carrier density in the GBL is rather small due to the 
compensation of the acceptor charges in the GBL and 
the donor and acceptor charges in the RIL (see, Sec.II). 
In this case, A 9 ~ eEj_dcBL , where E± = 47reE a/k 
is the transversal electric field created by the charged 
impurities and dcBL — 0.36 nm is the spacing between 
GLs in the GBL, so that A g ~ 4n e 2 T,AdGBL/n. Un¬ 
der the pumping conditions, the 2D electron and hole 
gases partially screen the electric field in the GBL struc¬ 
ture. As a result, one obtains A g ~ eEj_d e Jl L with 
d gbl < dcBL (35]. Taking into account the screening 
of the transversal electric field E± by both electron and 
hole components, we obtain the following estimate: A g ~ 
Ane 2 T,AdcBL /+ (8dGSi/as)], where as = «A 2 /me 2 
is the Bohr radius. Assuming as = 10 — 20 nm and 
E^ = (5 — 10) x 10 11 cm -2 , we obtain A 9 ~ 1.5 — 5 meV. 
These values are in line with those estimated in [35] and 
extracted from the experimental data (3, |45| . This esti¬ 


mate validates our model 

In the case of the clusterized charged impurities, the 
long range variations of the potential associated with the 
clusters lead to the variation of the band gap. The latter 
variations reinforce the electron and hole scattering on 
the clusters, but this effect is not really essential. 

In Eqs. (3) and (4), we have disregarded the indirect 
radiative processes associated with electron-hole scatter¬ 
ing in the GBL. The point is that the probability of such 
processes is proportional to [27r e 2 //t(g + QTi=’)] 2 E? 7 , where 
(in the degenerated two-dimensional electron and hole 
gases), the factors r] < 1 is the fraction of electrons and 
holes effectively participated in the scattering processes 
rj = T/£p at hio < T, so that Ery ~ E-r = 2mT/Trh 2 . 
The values of qTF and E?y are independent of the electron 
and hole densities and, hence, of the pumping conditions. 
This is because of the virtually constantdensity of state. 
At rooom temperatures, the value of E^ can be of the 
order of or less than E^ and much smaller than E ml- 
In the former case, one needs to replace the quantity D 
in equations (9) and (10) by D* = D/( 1 + 2Et/£a)- 

VIII. CONCLUSIONS 

We proposed to use the GBL-RIL heterostructure with 
the population inversion due to the electron and hole in¬ 
jection as the active region of interband THz lasers and 
demonstrated that the incorporation of the RIL enables 
a substantial reinforcement of the effect of negative THz 
conductivity and the laser THz gain. This is associated 
with the domination of the indirect interband transitions 
in the GBL over the indirect intraband transitions (re¬ 
sulting in the Drude absorption) when the carrier scat¬ 
tering on a long-range potential prevails. As shown, the 
latter can be realized due to the remote doping and en¬ 
hanced by the clusterization of impurities. 

Using a simplified device model, we calculated the THz 
conductivity of the GBL-RIL heterostructures as func¬ 
tion of the frequency for different structural parameters 
(dopant densities, spacer thickness, impurity cluster size, 
and the injected carrier momentum relaxation time). We 
found that the absolute value of the THz conductivity in 
the GBL-RIL heterostructures with a sufficiently highly 
doped RIL separated from the GBL by a spacer layer 
of properly chosen thickness can exceed in several times 
the pertinent value in the GBL-heterostructures without 
the RIL. Thus, the remotely doped GBL heterostructures 
can be of interest for applications as the active media in 
the THz lasers. 


Acknowledgments 

The authors are grateful to D. Svintsov and V. 
Vyurkov for useful discussions and information. The 
work was supported by the Japan Society for Promo¬ 
tion of Science (Grant-in-Aid for Specially Promoted 


Research # 23000008) and by the Russian Scientific 
Foundation (Project #14-29-00277). The works at UB 
and RPI were supported by the US Air Force award # 


FA9550-10-1-391 and by the US Army Research Labora¬ 
tory Cooperative Research Agreement, respectively. 


[1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K.S. 
Novoselov, and A. K. Geim, The electronic proper¬ 
ties of graphene, Rev. Mod. Phys. 81, 109 (2009), 
http://dx.doi.org/10.1103/RevModPhys.81.109 

[2] V. Ryzhii, M. Ryzhii, and T. Otsuji, Negative dynamic 
conductivity of graphene with optical pumping, J. Appl. 
Phys. 101, 083114-1-4 (2007), doi: 10.1063/1.2717566. 

[3] M. Ryzhii and V. Ryzhii, Ryzhii, Injection and popu¬ 
lation inversion in electrically induced p-n junction in 
graphene with split gates, Jpn. J. Appl. Phys. 46, L151 
(2007), doi: 10.1143/JJAP.46.L151. 

[4] V. Ryzhii, M. Ryzhii, V. Mitin, and T. Otsuji, 
Toward the creation of terahertz graphene injec¬ 
tion laser, J. Appl. Phys. 110, 094503-1-9 (2011), 
doi: 10.1063/1.3657853. 

[5] V. Ya. Aleshkin, A. A. Dubinov, and V. Ryzhii, Terahertz 
laser based on optically pumped graphene: model anf 
feasibility of realization JETP Lett. 89, 63(2009), doi: 
10.1134/S0021364009020039. 

[6] S. Boubanga-Tombet, S. Chan, A. Satou, T. Otsuji, 
and V. Ryzhii, Ultrafast carrier dynamics and tera¬ 
hertz emission in optically pumped graphene at room 
temperature, Phys. Rev. B 85, 035443-1-6 (2012), 
doi: 10.1103/PhysRevB.85.035443. 

[7] T.Otsuji, S. Boubanga-Tombet, A. Satou, H. Fukidome, 
M. Suemitsu, E. Sano, V. Popov, M.Ryzhii, and 
V.Ryzhii, Graphene-based devices in terahertz science 
and technology, J. Phys. D 45, 303001-1-9 (2012), 
doi: 10.1088/0022-3727/45/30/303001. 

[8] T. Li, L. Luo, M. Hupalo, J. Zhang, M. C. Tringides, 
J. Schmalian, and J. Wang, Femtosecond population in¬ 
version and stimulated emission of dense Dirac fermions 
in graphene, Phys. Rev. Lett. 108, 167401-1-4 (2012), 
http://dx.doi.org/10.1103/PhysRevLett.108.167401 

[9] Y. Takatsuka, K. Takahagi, E. Sano, V. Ryzhii, and 
T. Otsuji, Gain enhancement in graphene terahertz am¬ 
plifiers with resonant structures, J. Appl. Phys. 112, 
033103-1-4 (2012), http://dx.doi.org/10.1063/L4742998 

[10] T. Watanabe, T. Fukushima, Y. Yabe, S.A. Boubanga 
Tombet, A. Satou, A.A. Dubinov, V. Ya Aleshkin, V. 
Mitin, V. Ryzhii, and T. Otsuji, The gain enhancement 
effect of surface plasmon polaritons on terahertz stimu¬ 
lated emission in optically pumped monolayer graphene, 
New J. Phys. 15, 075003-1-11 (2013), doi:10.1088/1367- 
2630/15/7/075003. 

[11] T. Winzer, E. Marie, and A. Knorr, Microscopic mech¬ 
anism for transient population inversion and optical 
gain in graphene, Phys. Rev. B 87, 165413-1-5 (2013), 
http://dx.doi.org/10.1103/PhysRevB.87.165413 

[12] I. Gierz, J.C. Petersen, M. Mitrano, C. Cacho, I.C. Ed¬ 
mond Turcu, E. Springate, A. Stohr, A. Kohler, U. 
Starke, and A. Cavalleri,Snapshots of non-equilibrium 
Dirac carrier distributions in graphene, Nat. Mater. 12, 
1119 (2013), doi: 10.1038/nmat3757. 

[13] S. Kar, D.R. Mohapatra, E. Freysz, and A.K. 
Sood, Tuning photoinduced terahertz conductivity in 


monolayer graphene: Optical-pump terahertz-probe 
spectroscopy, Phys. Rev. B 90, 165420 (2014), 

http://dx.doi.org/10.1103/PhysRevB.90.165420 

[14] T. Li, L. Luo, M. Hupalo, J. Zhang, M.C. Tringides, 
J. Schmalian, and J. Wang, Femtosecond population in¬ 
version and stimulated emission of dense Dirac fermions 
in graphene, Phys. Rev. Lett. 108, 167401-1-4 (2012), 
http://dx.doi.org/10.1103/PhysRevLett.108.167401 

[15] R.R. Hartmann, J. Kono, and M.E. Portnoi, Terahertz 
science and technology of carbon nanomaterials, Nan¬ 
otechnology 25, 322001-1-16 (2014), doi:10.1088/0957- 
4484/25. 

[16] T. Otsuji, S. Boubanga-Tombet, A. Satou, M. Suemitsu, 
and V. Ryzhii, Spectroscopy study on ultrafast carrier 
dynamics and terahertz amplified stimulated emission in 
optically pumped graphene, J. Infrared Millimeter and 
Terahertz Waves, 33, 825 (2012), doi: 10.1007/sl0762- 
012-9908-8. 

[17] I. Gierz, M. Mitrano, J.C. Petersen, C. Cacho, I. C. E. 
Turcu, E. Springate, A. Stphr, A. Kphler, U. Starke, and 
A. Cavalleri, Population Inversion in Monolayer and Bi¬ 
layer Graphene, J. Phys.: Cond. Mat. 27, 164204-1-5 
(2015), doi: 10.1088/0953-8984/27/16/164204. 

[18] A. A. Dubinov, V. Ya. Aleshkin, M. Ryzhii, T. Otsuji, 
and V. Ryzhii, Terahertz laser with optically pumped 
graphene layers and Fabry-Perot resonator Appl. Phys. 
Express 2, 092301-1-3 (2009). 

[19] V. Ryzhii, M. Ryzhii, A. Satou, T. Otsuji, A. A. Du¬ 
binov,and V. Ya. Aleshkin, Feasibility of terahertz las¬ 
ing in optically pumped expitaxial multiple graphene 
layer structures, J. Appl. Phys. 106, 084507-1-6 (2009), 
doi: 10.1063/1.3247541. 

[20] V. Ryzhii, A. A. Dubinov, T. Otsuji, V. Mitin, and 
M. S. Shur, Terahertz lasers based on optically pumped 
multiple graphene structures with slot-line and dielec¬ 
tric waveguides, J. Appl. Phys. 107, 054505-1-3 (2010), 
doi:10.1063/1.3327212. 

[21] F. Rana, Graphene terahertz plasmon oscilla¬ 
tor, IEEE Trans. Nanotechnol. 7, 91 (2008), doi: 
10.1109/TNANO. 2007.910334 

[22] A. A. Dubinov, V. Ya. Aleshkin V. Mitin, T. Otsuji, 
and V. Ryzhii, Terahertz surface plasmons in optically 
pumped graphene structures, J. Phys. Cond. Mat. 23, 
145302-1-8 (2011), doi:1088/0953-8984/23/14/145302. 

[23] V. V. Popov, O. V. Polischuk, A. R. Davoyan, V. Ryzhii, 
T. Otsuji, and M. S. Shur, Plasmonic terahertz lasing 
in an array of graphene nanocavities, Phys. Rev. B 86, 
195437-1-6 (2012), doi:10.1103/PhysRevB.86.195437. 

[24] A. Tredicucci and M. S. Vitiello, Device concepts for 
graphene-based terahertz photonics, IEEE J. Select Top¬ 
ics in Quant. Electron. IEEE J. Sel Top. Quant. Electron. 
20, 8500109 (2014), doi: 10.U09/JSTQE.2013.2271692 

[25] H. C. Liu, C. Y. Song, A. J. SpringThorpe, and J. C. 
Cao, Terahertz quantum-well photodetector, Appl. Phys. 
Lett. 84, 4068 (2004), doi: 10.1063/1.1751620 . 

[26] B. S. Williams, Terahertz quantum-cascade lasers, Nat. 






9 


Photon. 1, 517 (2007), doi:10.1038/nphoton.2007.166. 

[27] S. Kumar, Recent progress in terahertz quantum cas¬ 
cade lasers, IEEE J. Sel. Top. Quantum. Electron. 17, 
38 (2011), doi: 10.1109/JSTQE.2010.2049735. 

[28] D. Svintsov, V. Ryzhii, and T. Otsuji, Negative 

dynamic Drude conductivity in pumped graphene, 
Appl. Phys. Express 7, 115101-1-4 (2014), doi: 

10.7567/APEX.7.115101. 

[29] D. Svintsov, T. Otsuji, V. Mitin, M. S. Shur, 

and V. Ryzhii, Negative terahertz conductivity in 
disordered graphene bilayers with population inver¬ 
sion, Appl. Phys. Lett. 106, 113501-1-4 (2015), 

http://dx.doi.Org/10.1063/l.4915314 

[30] M. Shur, Physics of Semiconductor Devices( Prentice- 
Hall, New Jersey, 1990). ISBN 0-13-666496-2 

[31] A. Shik, Quantum Wells: Physics and Electronics of Two- 
Dimensional Systems( World Scientific, Singapore, 1997), 
ISBN 981 02 3279 9. 

[32] T.Stauber, G. Gomez-Santos, F. Javier Garcia de 
Abajo, Extraordinary absorption of decorated un¬ 
doped graphene, Phys. Rev. Lett. 112, 077401 (2014), 
http://dx.doi.org/10.1103/PhysRevLett.112.077401 

[33] G. Konstantatos, M. Badioli, L. Gaudreau, J. Osmond, 
M. Bernechea, F. P. Garcia de Arquer, F. Gatti, F. H. 
L. Koppens, Hybrid graphene quantum dot phototransis¬ 
tors with ultrahigh gain, Nature Nanotechnol. 7, 363-368 
(2012), doi:10.1038/nnano.2012.60 

[34] E. McCann and V. Fal’ko, Landau-level degen¬ 

eracy and quantum hall effect in a graphite 
bilayer, Phys. Rev. Lett. 96, 0886805 (2006), 

http://dx. doi. org/10.1103/PhysRevLett.96.086805 

[35] E. McCann, D. S. L. Abergel, and V. I. Fal’ko,The 

low energy electronic band structure of bilayer. Eur. 
Phys. J. Special Topics 148, 91- 103 (2007),doi: 

10.1140/epjst/e2007-00229-l . 

[36] L. M. Zhang, Z. Q. Li, D. N. Basov, M. M. Fogler, 
Z. Hao, and M. C. Martin, Determination of the 
electronic structure of bilayer graphene from infrared 


spectroscopy Phys, Rev. B 78, 235408-1-5 (2008), 
http://dx.doi.org/10.1103/PhysRevB.78.235408 

[37] L. A. Falkovsky and S. S. Pershoguba, Optical far- 

infrared properties of a graphene monolayer and 
multilayer, Phys. Rev. B 76, 153410-1-4 (2007), 

http://dx.doi.org/10.1103/PhysRevB.76.153410 

[38] S. Das Sarma, E. H. Hwang, and E. Rossi, Theory of 
carrier transport in bilayer graphene Phys. Rev. B. 81, 
161407(R) (2010).doi: 10.1103/PhysRevB.81.161407 

[39] E. McCann, M. Koshino, The electronic properties of 
bilayer graphene, Rep. Prog. Phys. 76, 056503 (2013). 
doi: 10.1088/0034-4885/76/5/056503 

[40] H. Min, D. S.L. Abergel, E. H. Hwang, and S. Das 
Sarma, Optical and transport gaps in gated bilayer 
graphene, Phys. Rev. B 84, 041406-1-4(R) (2011), 
http://dx.doi.org/10.1103/PhysRevB.84.041406 

[41] F. T. Vasko and V. Ryzhii, Voltage and tem¬ 

perature dependencies of conductivity in gated 
graphene, Phys. Rev B 76, 233404 (2007). 

http://dx.doi.org/10.1103/PhysRevB.76.233404 

[42] C. H. Lewenkopf, E. R. Mucciolo, and A. H. Castro Neto, 
Numerical studies of conductivity and Fano factor in dis¬ 
ordered graphene Phys. Rev. B 77, 081410(R) (2008). 
http://dx.doi.org/10.1103/PhysRevB.77.081410 

[43] S. Adams, P. W. Brouwer, and S. Das Sarma, 

Crossover from quantum to Boltzmann transport 
in graphene Phys. Rev. B 79, 201404 (2009). 

http://dx.doi.org/10.1103/PhysRevB.79.201404 

[44] A. J. Samuels and J. D. Carey, Molecular Doping and 
Band-Gap Opening of Bilayer Graphene Alexander J. 
Samuels, and J. David Carey, ACS Nano, 2013, 7 (3), 
2790 (2013). doi: 10.1021/nn400340q 

[45] W. Zhang, Ch.-Te Lin, K.-ICu Liu, T. Tite, Ch.-Y. Su, 
Ch.-H. Chang, Yi-H. Lee, Ch.-W. Chu, K.-H. Wei, J.-L. 
Kuo, and L.-J. Li, Opening an electric band gap of bilayer 
graphene with molecular doping, ACS Nano 5 (9), 7517 
(2011).doi: 10.1021/nn202463g 




